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ABSTRACT 

From Zeeman Doppler images of the young, rapidly-rotating KO dwarf AB Doradiis, we 
have created a potential approximation to the observed radial magnetic field and have 
evolved it over 30 days due to the observed surface differential rotation, meridional 
flow and various diffusion rates. Assuming that the dark polar cap seen in Doppler 
images of this star is caused by the presence of a unipolar field, we have shown that 
the observed differential rotation will shear this field to produce the observed high- 
latitude band of unidirectional azimuthal field. By cross-correlating the evolved fields 
each day with the initial field we have followed the decay with time of the cross- 
correlatio n function. Over 30 days it decays by only 10%. This contrasts with the 
results of Barnes et al. (1998) who show that on this timescale the spot distribution 
of He699 is uncorrelated. We propose that this is due to the effects of flux emergence 
changing the spot distributions. 
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1 INTRODUCTION 



AB Doradus is a KO dwarf, orig inally discovered as a bright 
X-ray source by Pakull (1981)| an d subsequently observed 
by ROSAT ([Kurster et al. 1997| ) and BeppoSAX ( |Mag- 
gio et al. 200C ). Its photometric variability is believed to 



site polarities at intermediate to high latitude, which are 
approximately regularly spaced in longitude together with a 
unidirectional ring of azimuthal field at 70-80° indicatin g an 
underlying large-scale toroidal field (Donati et al. 1999). 



AB Dor exhibits flaring X-ray emission (Col 



this, cc mbined with its brightness (V ~ 6.8 — 7.0) and 199^ 

rapid rotation (P=0. 514) h ave made it an attractive can- 



be due to starspots (knders 1990i llnnis et al. 1988|) and licr Cameron et al. 1988 , gchmitt, Cutispoto fc Krautter 



didate for Doppler imaging (p^iirster, Schmitt fc Cutispoto 



199^ [Unruh, Collier Cameron fc Cutispoto 1995|,Ponati & 



1993 , ICollier Cameron fc Unruh 1994 ICollier C 



and there is indirect evidence from radio observations 
at 3,6,13 and 20 cm that the radio emissi on is highly direc - 
tive and suggests synchrotron radiation (Lim et al. 1994). 



Collier Cameron 1997 , Donati et al. 1999, Collier Cameron 
et al. 1999f)." From observations of the lithium line at 6708 A. 



The star is surrounded by a system of circumstellar promi- 
nences which can be observed as absorption transients in 
optically thick low-excitation lines e.g. H Balmer, Can and 
Mgll when the prominences cross the line of sight (Col 



lier Cameron fc Robinson 1989al Collier Cameron fc Robin-, 
son 1989byColli~ Cameron et al. 1990[ ). These prominences 



Rucinsl:i (1982 



suggested it was a post-T Tauri star. Ac- 



cording) to HIPPARCOS data it is 14.94±0.12 pc away 



Collier Cameron fc Foing (1997) inferred an age of ~2- 
3x10' year and its common-proper-motion companion, the 
M dwarf Rst 137b. 

AB Dor is of interest for a variety of reasons. The most 
important, for the purposes of this paper, is that Zeeman 
Doppler image s have been obtained in 3 consec utive years: 
1995 Dec 7-13(|Donati fc Collier Cameron 1997| ); 1996 Dec 
23-29( |Donati et al. 1999[ ); and 1998 Jan 10-15. These studies 
reveal that the radial field has at least 12 regions of oppo- 



E-mail: moira.jardine@st-and.ac.uk 



are trapped by the stellar magnetic field at, or beyond, the 
point of centrifugal balance. Their presence demonstrates 
that the corona is highly stru ctured even as far out as 3- 
5-R* (Jardine fc Ferreira 1996). 



Despite the rapid rotation, the differential rotation has 
been measured by cross-correlation of Zeeman Doppler im- 
ages secured a few days apart to be close to solar, with the 



equat or lapping the poles in ~110 ( i 
case) (Donati fc Collier Cameron 1997 



I 



120" in the solar 



onati et al. 1999) 



On longer timescales, however, Barnes et al. (1998) have 
shown that the spot distribution of the similar young rapid 
rotator He 699 becomes uncorrelated after 30 days. 
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The purpose this paper is to investigate the effect of dif- 
fusion and differential rotation on the evolution of AB Do- 
radiis's magnetic field. We aim to find out whether shearing 
at the edge of a unipolar cap can produce the observed ring 
of unidirectional azimuthal field. We also seek to determine 
the life times of surface magnetic features subject to diffusion 



and dif erential rotation. 



TEMPORAL EVOLUTION OF THE 
MAGNETIC FIELD 



filament channels on the Sun. It can also be used to study 
the field of AB Doradiis because we have high-resolution 
magnetic maps (2? 9 at the equator) and the differential ro- 
tation is similar to the solar value. The code takes the ob- 
served surface radial component of the field and calculates 
a potential field from this, and then evolves the calculated 
magnetic field due to the effects of differential rotation and 
diffiision^^^^^^^^^^^ 

Jardine et al. (1999) have demonstrated that for lati- 
tudes below about 60°the field is well-represented by a po- 
tential approximation. We anticipate that departures from 
a potential field caused by the shearing effect of the differ- 
ential rotation will appear at high latitudes near the edge of 
the unipol ar cap. By fitting the Stokes V profiles with both 
potential ( Hussain ^j^Jaxdm and 



non-potential (Hussain, van Ballegooijen & Jardine 2001) 
field models, it is possible to show that any currents are 
concentrated close to the pole. 



2.1 The scalar magnetic potential, i/) 

If we assume that the field is potential, then we can write 
B in terms of a fiux function ip, with 

B = ~ViP, 

which in spherical co-ordinates gives 

dip 



Br 



Bo 



and 



-B<4 



dr ' 

1 9i/> 
'r de 

1 dtp 



r sin 9 d(f> 

Here ip satisfies Laplace's equation V^V ~ which can be 
expressed as 

1 9V 



1 d , 2dip 



1 



dr d: 



■sine 89 



sin^ 9 dcj)'^ 



= 0.(1) 



A separable solution for i/) can be found 

JV I 
i — 1 m— — Z 

where Pim are the associated Legendre functions and 

tpim{r) = ai„y + bimr~^''^^\ 



We chose to truncate the series at A'' = 63, corresponding to 
the maximum resolution of the reconstructed field images. 
We clearly need two boundary conditions to determine 
We chose to specify as one boundary condition that at some 
distance from the star (the source surface, r, ~ 5.1-R, ), the 
field is radial and so Bg{rs) = B^{rs) — (Schatten, Wilcox 



Ness 1969). This mimics the stellar wind. 



Since most stellar prominences form at around the 
corotation radius (2.7R*), we know that a significant frac- 
tion of the field is closed at that radius. Hence we choose 
rs — 5.1ii*. We then have 



We are using a code originally d eveloped by van Ballegooi- 
jen, Cartledge & Priest (1998) to study the formation of equivalent to 



ai; 



i-1 I J. 
^r^ + bimrs 



= 0, 



0. 



As a second boundary condition we impose the radial field 
at the surface to be the observed radial field. We can then 
express the magnetic field in terms of the two-dimensional 
Fourier coefficients Bim, where 

Bira{R,)^2'K [ B^{9)Pi^{9) sin 9d9 



Blm{R,) = -laimRi ' + (/ + l)hl^R^'-^. 

The function Bm {&) is derived from a fast Fourier transform 
performed latitude-by-latitude on the observed radial field 
Br{R,,9,<l,). 

Once the field is evolved due to diffusion and differen- 
tial rotation, it is not necessarily potential, although it can 
still be expressed as a sum of spherical harmonics. The field 
components are then expressed in terms of the functions 



{9)e 



l — l m,— — l 



and 



where 



l — l m,— — l 



J 



-ism9B,)- — 



sin 9 

is the radial component of the current and 
1 



A = 



d , . an \ , dB^ 



is the 2-dimensional divergence. 

2.2 Evolving the field using the induction 
equation 

From three of Maxwell's equations 

VxB = ^ij, 

V.B = 0, 
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Figure 1. This shows how the magnetic energy varies as the field evolves for 1998. The left hand diagram is for the case when 
r;=250km^s~^ , with differential rotation (solid line) and without (dotted line), and the right hand diagram compares the case where 
there is no diffusion, i.e. r) = (solid line) and where r7=250km^s~^ (dotted line). 




Figure 2. Cross-correlation of calculated radial field with the observed radial field for 1995, latitude 30°(left) and 60°(right). The solid 
line represents rj = 250km'^s~^, the dotted line rj = 350km^s~^, the dashed line rj = 450km^s~^ and the dot-dash line r} = 550km^s~^. 
We should not expect a cross-correlation of exactly 1 at day since we are correlating an observed field with a calculated one. 



and Ohm's Law 



j = cr(E + vxB) 



where a is tho conductivity wo got tho induction oquation 



^ = Vx(vxB)-VxE'. 
Here E' is given by 



(2) 



V 



r sin 9 



-(B.sm^)- — 



— ; T] — 

-tifl = — : — rr^TT- 



Here A = -I — ^ is the latitude, Ao gives the latitude above 
which the meridional flow is zero, Ao = 75° and uq = 11m 
s~^ which is close to the predicted value ( Kitchatinov & 
Riidigor 1999). The values of Bjm, ^im and Jim are evolved 



using the induction equation according to the meridional 
flow, the observed differential rotation and using various 
values of the magnetic diffusion, ranging from 250 to 550 
km^s~^ (c/. the solar value of 450 km^s~^). 
The differential rotation is of the form 

n(6l) = 12.2434 - 0.0564 cos^ 9 rad d (4) 



wher e Q is the rotation rate (Donati & Collier Cameron 



1997). We have assumed that r/ is uniform across the surface, 
although there is evidence t hat this may not be the case for 
the Sun (|Berger et al. 199^). 



1] dBr 



with ri =l/^ia being the mag netic diffusivity (van Ballegooi 
jen, Ca'tledge & Priest 1995). This assumes that there is no 
radial transport of the magnetic field, and that the merid- 
ional flow is poleward and the same as the solar value given 
by 



-uosin(5^) ff |A| < Ao 
I otherwise. 



(3) 



2.3 Calculating the evolved field 

The third stage is to take the evolved coefficients Bim.,Aim 
and Ji,n and the associated Legendre functions Pim and 
calculate the three components of the magnetic field- Br 
(radial), B^ (azimuthal) and Be (meridional)- from them. 
These will be given by 

JV I 
l — l m— — l 
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Figure 3. This diagram shows the azimuthal and radial components of the field in 1996. The top row gives the observed radial (left) 
and azimuthal (right) fields for Dec 23/25. The middle row gives the observed radial (left) and azimuthal (right) fields for Dec 29. The 
bottom row gives the calculated radial (left) and azimuthal (right) fields for Dec 29, assuming that )7=450km^s~^ 



Be{r,e, 
where 



l — l m—~l 
N I 



. , ^ dPim , . , mPim 

Alm(r) h tJlm . „ 

dO sm 



imPim _ dPim 

^) . /I -Jim 

inO dO 



imip 



l — \ m— — l 



3 RESULTS 

As an initial consistency check we computed the evolution 
of the magnetic energy in the field at the surface, calculat- 
ing the ratio of magnetic energy in the evolved case to the 
original. 

A magnetic field has energy _B^/2/i per unit volume, so 
the total energy is 



W 



dV 



volume ' 

Fig |l| shows this for the f998 field. The left-hand panel 
shows the degree of diffusive decay of the field energy over 



30 days. The right-hand panel shows the increase in energy 
that would occur in the absence of diffusion, due to the 
winding-up of the field by the differential rotation. 



3.1 Long-term evolution of the surface field 

We took the observed field for the 1995 and 1998 observa- 
tions and evolved it over 30 days according to the induction 
equation (4). This allowed us to study the variation with 
time of the cross-correlation of the radial component ob- 
served on the first night with that calculated for subsequent 
nights (Fig We chose two latitudes: 30 and 60° north. 
The results for 1998 are qualitatively similar. In all cases, 
the cross-correlation function decays by approximately 10% 
over 30 days. Although choosing a higher value for the diffu- 
sivity does cause a more rapid decay of the field and hence 
a faster decay of the cross- correlation function, it is still not 
en ough to explain the complete lack of correlation found 



by Barnes et al. (1998) for He699. It appears that for AB 
Doradtis, if diffusion and differential rotation were the only 
processes causing the field to evolve, that even after one 
month there should still be a good correlation. 
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Figure 4. The peak amplitude of the cross-correlation function between the observed and calculated radial (left) and azimuthal (right) 
fields for latitudes between and 80°north for 1996 Dec 29. The solid line represents r] = 250km^s~^, the dotted line r] = 350km^s~^, 
the dashed line r] = 450km^s^^ and the dot-dash line r] = 550km^s^^. 





Figure 5. The peak amplitude of the cross-correlation function between the observed and calculated radial (left) and azimuthal (right) 
fields for latitudes between and 80°north for 1996 Dec 29 with rj = 450km'^s~^. The solid line is with the differential rotation, and the 
dotted line is without the differential rotation. 



3.2 Short-term evolution 

For comparison, we can look at the field evolution over a 
much shorter period of time. This has the advantage that 
we can compare our results with the observed evolution of 
the field during a single observing run. Here we use data 
from the 1996 run, in which we secured two sets of mag- 
netic maps, separated by 5 nights (Fig. Our aim is to 
determine whether the observed magnetic elements retain 
their identities over a period of 5 nights in the presence of 
diffusion and differential rotation. 

We began by investigating the effect of varying the dif- 
fusivity. We evolved the 1996 Dec 23/25 field forward in time 
using various values of the diffusion coefficient. In each case 
we cross-correlated the resulting radial field map with the 
observed radial field for Dec 29 (Fig We ensured that the 
cross-correlation only invo lved those longitud es that were 
well-observed, viz. 18-180° (Jardine et al. 199£). The cross- 
correlation function was computed for each of a set of lati- 
tudes between and 80° in the visible hemisphere, and the 
amplitude of the strongest peak in the ccf at each latitude 
is plotted in Fig. ^. It can be seen that the five-day span of 
these observations is too short for differences in the value 
chosen for diffusivity to have much effect. 

We also considered the case where the field was allowed 
to evolve under the influence of diffusion only, i.e. with the 



differential rotation switched off. We found that over the 5- 
day span of the 1996 December observations, the influence 
of the differential rotation was negligible (Fig. 

From these we see that altering the value of the diffu- 
sion and removing the differential rotation have little effect 
on the cross- correlation over 5 nights. We would need obser- 
vations over a longer timescale, say a month, to be able to 
look for meaningful results. 



4 THE HIGH-LATITUDE AZIMUTHAL FIELD 



Jardine et al. (1999) demonstrated that the high-latitude az- 
imuthal band of held was not reproduced by modelling the 
fleld as potential. Here we investigate whether the differen- 
tial rotation could produce this band by taking the 1998 
radial fleld map and adding in a cap of unipolar radial fleld 
extending from the pole to latitude 80° , well within the dark 
polar spot seen in the stellar surface-brightness map. An 
identical cap of opposite polarity was added to the unob- 
servable hemisphere to conserve flux. We emphasize that the 
polarity and strength of the dark polar cap cannot be deter- 
mined directly from the observations, since the low surface 
brightness and strong foreshortening suppress the Zeeman 
signal from this part of the star. For each of a range of plau- 
sible polar field strengths and polarities we evolved the field 
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forward in time for 5 days, then compared the mean value 
of at each latitude with the observed value. 

We see from Fig. ^ that a high-latitude azimuthal band 
of negative polarity is produced when we shear the image 
at the observed differential rotation rate with a polar cap 
having Br < Q and vice versa. The results confirm that 
the observed differential rotation is capable of producing 
a high-latitude negative azimuthal band, as is seen in the 
observations from all three observing seasons (Fig.|). This 
predicts that the magnetic polarity of the dark polar region 
was predominantly negative in all three seasons. 

The latitude of the maximum in B^ occurs at the edge 
of the imposed unipolar cap in Fig. ^ The azimuthal field 
is localised here because most field lines originating in the 
north polar cap are connected to high northern latitudes just 
outside the cap. The direction of these field lines depends 
strongly upon the distribution of radial field at mid to high 
latitudes. Changing the strength of the polar cap has little 
effect on field lines at low latitude. This corresponds to our 
result that ^ B^ varies little at low latitude as the strength 
of the cap is altered (as seen in Fig |^ . 

The observed azimuthal band plotted in Fig. |^ is more 
diffuse than the model, having a broader peak between lati- 
tudes 65° and 80° in all three seasons' data. This corresponds 
roughly to the edge of the dark polar region, which extends 
to latitude 70°or so. The breadth of the peak suggests a 
more gradual fall-off in the polar field than we imposed on 
the model. The apparent decrease in field strength at lati- 
tudes above 70° can probably be ascribed to suppression of 
the Zeeman signal within the dark polar region. 

The strength of the azimuthal field after only 5 nights 
is substantially less than that observed. This is not surpris- 
ing, since the timescale on which shear can generate an az- 
imuthal field with a strength comparable to the radial and 
meridional field will be of order the equator-pole lap time 
of 110 days. The diffusion time for length scales comparable 
to the size of the individual magnetic regions in the images 
is also of this order. Once an equilibrium is established be- 
tween diffusion and differential rotation, we would expect 
the azimuthal field strength to be an order of magnitude 
greater than that produced after 5 days, in agreement with 
the observations. On a 100-day timescale, however, we ex- 
pect the picture to be complicated further by the emergence 
of new flux, making a direct comparison with the observa- 
tions problematic. 



5 CONCLUSIONS 

We have modelled the evolution of the magnetic field of 
AB Doradiis due to the effects of differential rotation and 
diffusion. We use as a starting point the Zeeman Doppler 
images obtained on three consecutive years and assume that 
the field is initially potential, but evolves away from this 
state as a function of time. 

Over a timescale of 20 to 30 days we have deter- 
mined, as a function of time, the cross-correlation of our 
model radial magnetic field with the observed radial com- 
ponent on the first night. We find that over one month, 



vulje of the evolved B, for 1998 




Figure 6. Azimuthally-averaged values of as a function of 
latitude computed from the 1998 data, evolved for 5 nights with 
ri = 250km^s~^ and with different values of the polar field. The 
solid line is the case where the polar field is set to OG, the dotted 
line where it is set to +1000G, the dashed line where it is set 
to -lOOOG, the dash-dot line where it is set to -I-4000G and the 
long-dashed line where it is set to -4000G. 
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Figure 7. Azimuthally-averaged values of as a function of 
latitude in the observed field maps for 1995 (solid line), 1996 Dec 
23/25 (dotted line) and 1998 (dashed line). 



this timescale gives much more rapid decrease of the cross- 
correlation function. This result suggests that the evolution 
of AB Doradhs's surface magnetic field is not governed solely 
by diffusion and differential rotation. We conclude that these 
results are more likely to be due to the effects of flux emer- 
gence changing the spot distribution than the effects of dif- 
fusion or differential rotation. 

This result is independent of the assumed degree of fleld 
diffusion. We have compared the effects of values of r) rang- 
ing from 250 to 450 km'^s"^ and found the results to be 
qualitatively the same. The presence of some diffusion is of 
course necessary (and we have confirmed that the magnetic 
energy grows monotonically with time in the absence of dif- 
fusion). The exact value of r] seems however to have little 
effect on the five-day timescale of a typical observing run. 
Since diffusion has little effect on the flux distribution, the 
differential rotation acts simply to advect the field. Con- 
sequently, although at each latitude the peak of th e cross- 



correlation function may be at a different longitude ( Donati 



the cro ss- correlation function decays by about 10%. Ob- 



servations of He699 by Barnes et al. (1998) show however 
that cross-correlating the observed spot distributions over 



& Collier Cameron 1997), its actual value is virtually un- 
changed by the effects of differential rotation. 

We have also compared the radial and azimuthal mag- 
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netic fields generated by our model over 5 nights with those 
obtained from Zeeman Doppler images on 1996 Dec 23/25 
and 29, and found the agreement to be excellent. The evo- 
lution of the azimuthal field is of particular interest with 
regard to the band of high latitude unidirectional azimuthal 
field seen in the Zeeman-Doppler images. We have investi- 
gated whether the shearing effect of the differential rotation 
is sufficient to generate this band of field. The polaxity of 
this band depends upon the sign of the radial field in the 
polar cap, whereas its strength depends on the competition 
between shear and diffusion. Since the diffusion timescale 
for resolvable features is comparable to the winding time, 
wc expect the azimuthal field to attain a strength compa- 
rable to the radial and meridional field near the boundary 
of the polar cap, as is indeed observed. Our results suggest 
that the differential rotation could play a major part in the 
creation and preservation of a high latitude azimuthal band. 
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